Vortices in fermion droplets with repulsive dipole-dipole interactions 
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Vortices are found in a fermion system with repulsive dipole-dipole interactions, trapped by a 
rotating quasi-two-dimensional harmonic oscillator potential. Such systems have much in common 
with electrons in quantum dots, where rotation is induced via an external magnetic field. In 
contrast to the Coulomb interactions between electrons, the (externally tunable) anisotropy of the 
dipole-dipole interaction breaks the rotational symmetry of the Hamiltonian. This may cause the 
otherwise rotationally symmetric exact wavefunction to reveal its internal structure more directly. 
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I. INTRODUCTION 

Rotating quantum fluids have been studied for a long 
time, with the interest being spurred by their many fas- 
cinating properties. A prominent example are superflu- 
ids such as liquid helium [lj. When set rotating in a 
bucket, for sufficiently large rotation, the quantum liq- 
uid becomes penetrated by quantized vortices, forming 
the well-known Abrikosov lattice. 

After the early prediction of Bose and Einstein in the 
1920's [2 J it was not until 1995 that the condensation of 
a gas of bosonic atoms into a single coherent quantum 
state could be achieved [3J. Stirring the condensate with 
lasers or rotating the trap that confines the dilute atom 
gas, vortices similar to those in 4 He [4j [5] were observed 
[61 [7] (see also the review by Fetter [8]). Similar states 
have been realized also for trapped fermionic atoms with 
attractive interactions where pairing or molecule forma- 
tion can occur (see the reviews [9) HO])- 

Apart from their presence in systems with bosons, vor- 
tices have also been predicted to occur in fermion systems 
with purely repulsive interactions, such as quantum dots 
- small man-made electronic systems that can be created 
in a semiconductor heterostructure [33], where the rota- 
tion can be induced with an external magnetic field. For 
small numbers of confined electrons, it was shown that 
these fermionic quantum Hall droplets form vortices in 
a very similar way than repulsive bosons set rotating in 
the trap [TT1IT3] . This analogy may, in fact, question the 
commonly accepted view that vortices and vortex arrays 
may be taken as a criterion of superfluid properties [14J. 
Experimentally however, vortices in quantum dots are 
difficult to detect: As the electrons are inside a semicon- 
ductor crystal, probing their properties must typically be 
done by indirect methods. Examples are electron trans- 
port or magnetization measurements [TU [15]. These are, 
however, strongly hindered by the restricted resolution in 
the conductance spectra, as well as unavoidable sample 
imperfections. 

Here, ultra-cold atomic gases may be the better choice, 
as they typically are very clean, and remarkable charac- 
terization techniques have been demonstrated. For ex- 
ample, it is possible to directly image the atomic cloud 



after expansion. Atomic quantum gases usually confine 
millions of atoms, bringing system sizes close to the ther- 
modynamic limit. Serwane et at. [16] however showed 
recently that the confinement of about a dozen of cold 
fermionic atoms can be reached experimentally. 

The electrons in a quantum dot interact via long- 
range repulsive Coulomb forces. An alternative is pro- 
vided by atoms or molecules with (either electric or mag- 
netic) dipole-dipole interactions providing a long-range 
coupling between the particles. (There has recently been 
much interest in such dipolar systems, see for example 
Refs. [T7H22] , and [23 ) [24 ] for reviews). In contrast, the 
very short-ranged van der Waals interaction has a limited 
effect on spin-polarized fermions due to the Pauli exclu- 
sion principle, and might not give a similar response to 
rotation as long-range interactions. 

In this paper, we show that vortices and vortex clusters 
occur in a (spin-polarized) fermion droplet with strictly 
repulsive dipole-dipole interactions, in much analogy to 
vortices in a quantum dot. While in quantum dots the 
vortex lattice still awaits experimental detection, we sug- 
gest that dipolar fermionic atoms with repulsive inter- 
actions may indeed show a similar vortex lattice than 
analogous systems with bosons. 



II. MODEL 

As a simple model (that well resembles that of elec- 
trons in a quantum dot) we consider a few (spin- 
polarized) fermions confined in a rotationally symmet- 
ric two-dimensional harmonic oscillator in the x?/-plane, 
with oscillator frequency ujq. We assume that the par- 
ticles are confined in the z-direction with a tight har- 
monic oscillator, so that the total trapping potential is 
VtTup(x,y,z) = \mujl(x 2 + y 2 ) + \muo 2 z z 2 . The oscil- 
lator length l z is here set to be 1/100 of that in the 
plane, so that the particles occupy the lowest orbital 
in this direction. (In the following, we work in dimen- 
sionless os cillator u nits, e.g. lengths are given in units 
of Iq = yfh/ (mcjo)? energies in Hloq and frequencies in 
ljq.) The rotation is induced by adding a term —QL to 
the Hamiltonian, where Q is the rotational (angular) fre- 
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quency of the rotating trap, and L is the (z-projection) 
of the angular momentum operator. 

For the interaction between the fermions we assume 
that their dipole moments are aligned by an external 
field. This makes it possible to control the effective 
anisotropy of the interaction, breaking the rotational 
symmetry in the system. A few other theoretical studies 
considered very similar setups [25-27J. Here, however, 
the focus is on the vortex structure. 

For two particles with dipole moments, a general ex- 
pression for their interaction energy is given in e.g. Ref. 
[23]. In the present system, we assume the dipole mo- 
ments to be aligned to an axis lying in the xz-plane, such 
that this axis forms an angle with the x?/-plane where 
the particles are effectively confined by the harmonic 
trapping. For two point dipoles confined to the x7/-plane, 
the interaction is attractive for 6 < arccos ^ « 54.7°. 
To avoid the possibility that the system collapses in the 
case of attractive interactions, we restrict the tilting to 
the interval 54.7° < B < 90°, where the interaction is re- 
pulsive. The coordinate system used here is co-rotating 
with the trap (at frequency f2), meaning that the dipole 
axis is rotating in the laboratory frame. The in-plane in- 
teraction between the particles is obtained by analytically 
integrating their motion along the z-direction, where the 
particles are assumed to be in the lowest oscillator state. 
This results in the interaction potential 

D 2 e^ 2 

^2D(r, </>) =^= {(2 + 20#o(£/2) - 2^1 (e/2) 

+ cos 2 6[-(3 + 2£)K (£/2) + (1 + 20 
x Ki(£/2)] +2cos 2 Bcos 2 0[-^ o (^/2) 
+ (£-l)ifi(£/2)]}, 

(1) 

where K$ and K\ are irregular modified Bessel func- 
tions, and £ = r 2 /(2(I z /Iq) 2 ). The prefactor D is here 
dimensionless, for an electrical dipole it corresponds to 
D = J^- — -^^= where d is the dipole moment, and simi- 

V47reo ^V^o 

larly for a magnetic dipole [31 . The expression in Eq. ([!]) 
(or special cases of it) is given e.g. in Refs. (26[ I28l43l[ . 
Only for = 90° does the interaction have rotational 
symmetry. In this limit, it is qualitatively similar to the 
electrostatic Coulomb interaction, however, being more 
short-ranged. For other angles, the two-body term is 
spatially anisotropic (cf Fig. 1 in Ref. |31j). 

The ratio l z /lo enters above, and as it is typically much 
smaller than one, at first sight it gives a very large inter- 
action strength. However, this ratio also enters elsewhere 
in the expression so that its final behavior is more subtle. 
It turns out that for long ranges, Eq. ([!]) is practically 
independent of I z /Iq and proportional to D 2 /r 3 . (Long 

range here corresponds to r ^> l z /h-) At short range 
n 2 

the coefficient ^ /i )z does have an impact for the value 
of the interaction potential, but in this study we only 
consider spin-polarized fermions, for which the Pauli ex- 
clusion principle eliminates the effect of the interaction 
at short ranges. 



The many-body Hamiltonian 
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is then diagonalized numerically, applying the configura- 
tion interaction method (also called 'exact' diagonaliza- 
tion). The commonly used lowest Landau level approx- 
imation [14] allows us to restrict the Hilbert space to a 
basis of Slater determinants constructed with the single- 
particle orbitals ip n =o,m>o(r, <p) = ^_L_ r m e zm ^e~ r I 2 of 
the two-dimensional harmonic oscillator. This is a rea- 
sonable approximation in the limit when the rotational 
frequency ft is close to, but smaller than, the trapping 
frequency c^o, and the interaction is weak. We set the 
interaction strength D = 0.1 to assure that the inter- 
action energy is smaller than the Landau level spacing, 
which is of order Hujo. In the calculations, sufficiently 
large m- values are included so that the resulting energies 
are unaffected by this cutoff. 

III. RESULTS 

Characteristic for the response of the system to a 
monotoneously increasing trap rotation Q is the change 
in angular momentum L(f2), obtained by minimizing the 
total energy in the rotating frame of reference as a func- 
tion of ft. For a supernuid, after the onset of trap ro- 
tation, the system first remains at rest, L = 0, until a 
critical frequency is reached beyond which vortices begin 
to penetrate the cloud. As Q further increases, the an- 
gular momentum changes discontinuously to higher val- 
ues as additional vortices appear. For a rotating Bose- 
Einstein condensate, this sequential appearance of vor- 
tices was predicted by applying the mean-field Gross- 
Pitaevskii approach [4, 5j, and later also observed ex- 
perimentally [6j[7]- 

For electrons in quantum dots, the rotation may be 
induced by an applied magnetic field. In the limit of slow 
rotation, vortices and vortex clusters occur even though 
there is no condensation of bosons present. (For a more 
detailed discussion, see e.g. Refs. [T2J E]). 

Turning now to fermions with repulsive dipole-dipole 
interactions, for small particle numbers and relatively 
slow rotation considered here, we find that indeed, the 
system shows a similar vortex pattern. 

Solving the many-body Hamiltonian by direct diago- 
nalization, the numerical effort grows considerably with 
particle number. Because of the anisotropic dipole-dipole 
interaction, angular momentum is not a good quantum 
number, which poses an additional limitation to the sys- 
tem sizes that are numerically accessible. For this reason 
we here limit the study to six particles confined in the 
harmonic trap, making use of the Lowest Landau level 
appropriate for weak interactions. 

In the case of non-interacting bosons, the many- 
particle ground state is the permanent |7V0000 . . . ) with 



3 



macroscopic occupancy of the m = orbital in the LLL 
basis. For spinless fermions, however, the Pauli princi- 
ple demands single occupancy of the orbit als, leading to 
the so-called "maximum density droplet" 1 111111000 . . . ). 
This state represents the finite-size analog to the Laugh- 
lin state for the integer quantum Hall effect, and is the 
fermionic equivalent to the condensate permanent. The 
lowest possible angular momentum of this Slater deter- 
minant is Lmdd = N(N — l)/2. 

The response of the system with TV = 6 dipolar 
fermions to the increasing trap rotation Q is shown in 
Fig. [T] for different values of the dipolar tilt angle where 
the interactions are in the repulsive regime. (Note that 
since L is not a good quantum number, we instead cal- 
culate the corresponding expectation value that can be 
non-integer.) Qualitatively, the resulting picture is very 
similar as for trapped bosons, or electrons in quantum 
dots. 



A. Isotropic interaction 

Let us first consider the case of aligned dipoles with 
= 90.0°. At moderate rotation, the cloud remains at 
rest relative to the angular momentum of the MDD (here, 
Lmdd = 15), until a critical frequency is reached, where 
the first step in L(Q) occurs, and the system aquires an- 
gular momentum beyond the MDD value. At this fre- 
quency a vortex is formed at the center of the cloud 
(see discussion below). With increasing rotation, fur- 
ther steps occur that are associated with the entry of 
additional vortices. The tilt angle of the dipole axis is 
found to have little effect on L(Q): The jumps are mainly 
shifted to higher values of Q as the angle is lowered. This 
can be understood by noting that the tilt of the dipole 
angle not only makes the interactions anisotropic, but 
also effectively lowers their strength. 

The single-particle density distribution, however, de- 
fined as p(r) = (^|^ t (r)^(r)|^), can be strongly af- 
fected by the tilt, as shown in Fig. [2j If the dipoles are 
aligned perpendicular to the plane of motion (0 = 90°), 
the Hamiltonian has rotational symmetry. Hence, the 
density of the associated eigenstates should also be az- 
imuthally symmetric. This makes it impossible to inden- 
tify vortex clusters beyond the simple unit vortex (that 
has azimuthal symmetry). Internally, for the two- vortex 
solution, the two single vortices appear as a pair with 
two-fold symmetry, while the three-vortex solution lo- 
cates the vortices at the corners of an equilateral trian- 
gle. However, if the Hamiltonian has rotational sym- 
metry these structures cannot be visible in the single- 
particle density. (The pair-correlated density, though, 
can reveal them, as discussed later in the text.) This is 
in contrast with wavefunctions obtained from mean-field 
approximations such as the Hartree-Fock or the Gross- 
Pitaevskii equations, where broken-symmetry solutions 
are possible. 
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Figure 1. (Color online) Expectation value of the angular 
momentum L of the ground state, as a function of the rota- 
tional frequency Q, for three different tilt angles of the dipoles. 
The system consists of N — 6 (spin-polarized) fermions with 
dipole-dipole interactions. As Q increases, the angular mo- 
mentum changes discontinously, and vortices penetrate the 
quantum system. The states marked with circles are the ones 
which are shown in the following figures. 



B. Anisotropic interaction 

When the dipole axis is tilted, the rotational symme- 
try is broken. Though some states are not significantly 
affected by the tilt - e.g. the state with a single vortex at 
the center - with increasing ^, the state with (L) « 25 
now shows two clear off-center minima in its density (see 
Fig. [2|. The overlap of this state at = 55.6° with that 
at = 90° is 93.1%, showing that the internal structure 
of the quantum state is largely unchanged despite the 
seemingly different densities. 

The density does not drop to zero at the vortex cores, 
but this typically does not happen in a finite-size system. 
For example, rotating bosons forming a single vortex ex- 
hibits zero density in the mean-field solution, but has a 
nonzero density at the vortex core in the exact analytical 
solution for a specific finite number of particles [34H36] . 



C. Currents 

While local minima in the density may indicate vor- 
tices, further support for their existence can be obtained 
by examining the probability current, as shown in Fig. [3] 
The current is here given as the expectation value of the 
operator (cf. equations 3-5 in Ref. [32] ) 

N 

j(r) = J2 y [<K r - r Vi + Vi$(r -Ti)}- (Sle z x r)p(r) 

2 = 1 

(3) 
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(L) = 90.0° = 72.8° = 55.6° 
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Figure 2. (Color online) Probability densities for the states 
marked with circles in Fig. [I] The listed angular momentum 
values are approximate. Each subplot, here and in subsequent 
figures, is plotted in the intervals — 4 < x < 4 (horizontally) 
and —4 < y < 4 (vertically). All plots are normalized to have 
the same peak height, and plotted according to the shown 
colorbar. 

where p(r) — J^Li — an d &z is the unit vector 
of the z-direction. The last term in the expression origi- 
nates in the coordinate transformation we use here - the 
current (as all other calculated quantities) is given in the 
co-rotating reference frame. In principle one could obtain 
the velocity v by dividing the current with the density, 
and then the vorticity as V x v. However, in the regions 
where the density is very small this could yield numer- 
ically unstable results, and we here restrict ourselves to 
analyzing the more well-defined current. 

For the two- vortex state discussed earlier, Fig. [3] shows 
the corresponding current which can be seen to circulate 
around the density minima, supporting the interpreta- 
tion of these minima as vortices. We note that in Fig. [3] 
the inner and outer regions of the system appear to be 
rotating in opposite directions - this is however an effect 
of the co-rotating reference frame used here. 



D. Pair-correlated densities 

In Fig. [4] we show pair-correlated densities. Comple- 
menting the single-particle density, this quantity can give 
additional information about the internal structure of 
the state. The pair-correlated density is here defined as 
p(r,r') = (^|^ t (r)^ t (r / )^(r / )^(r)|^), giving the prob- 




Figure 3. (Color online) Probability current (as defined in the 
text) for the case (L) « 25, when G = 55.6°. The inset shows 
the case (L) « 21, as comparison. The black arrows show the 
current, on top of the single-particle density. In both cases 
shown, the current can be seen to loop around the density 
minima. As an effect of the calculations being performed in 
the co-rotating reference frame, the current in the outer parts 
of the particle cloud appear to rotate opposite to the current 
around the vortices. The same plotting conventions as in 
Fig. [2] are used here. 



ability density of simultaneously finding two particles at 
positions r and r' . In Fig. [5] we show pair-correlated 
densities for the state at = 55.6° with a reference par- 
ticle is placed at different positions r' - showing that the 
overall structure of the density is a stable configuration 
with two minima, although somewhat obscured here by 
the so-called exchange hole around the position of the 
reference particle. (This effect is due to the Pauli exclu- 
sion principle, which implies that the probability to find 
any other particle must vanish near the position of the 
reference particle.) 

For larger rotation, as ft is increased further, the N = 
6 system appears to have too few particles to support 
more vortices, and instead it undergoes a transition to a 
state where the particles are highly localized, as shown 
m Fig. |4| This is again very similar to the behavior of 
electrons in a strong magnetic field. Though not seen for 
the parameter regimes considered here, the anisotropic 
interaction can be expected to affect the geometries of 
such states (see Ref. |27J). 
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Figure 4. (Color online) Pair-correlated densities for the 
states shown in the right column of Fig. [2] for which the dipole 
tilt angle is G = 55.6°. The white cross marks the position 
of the reference particle. The same plotting conventions as in 
Fig. [2] are used here. The positions of the reference particle 
are here all on the x-axis, and from left to right x =1.5, 1.8, 
2.3 and 2.4. 
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Figure 5. (Color online) Pair-correlated densities for the state 
with (L) « 25, for G = 55.6°, with the reference particle 
placed at different positions (white crosses). Apart from the 
necessary fermionic exchange hole, the general structure of the 
state is not sensitive to the chosen position. The positions 
of the reference particle, starting from top left panel, were 
(x,y) =(2.4,0), (2.2,0.6), (1.9,1.2), (1.5,1.5), (0.9,1.7) and 
(0, 1.7). The same plotting conventions as in Fig. [2] are used 
here. 



IV. SUMMARY 

To summarize, we have shown that a quasi-2D trap 
with fermionic particles with dipole-dipole interactions 



show a response to induced rotation that is in many ways 
similar to that of confined electrons in a magnetic field 
|14j . This is perhaps naturally expected given that both 
the appearing interactions in both cases are repulsive and 
long-ranged (as also discussed earlier in Refs. (25) [26]) . 
Here, however, we focused on the occurrence of vortices 
in repulsive fermion systems, that may show very similar 
properties than their superfluid cousins with bosons or 
attractive fermions. We found that a system of six dipo- 
lar fermions in a harmonic trap shows the characteristic 
step-wise increase of angular momentum as a function 
of the trap rotation. While the cloud first remains at 
rest, beyond a critical frequency the angular momentum 
shows a discontinuous jump by N units, and the first vor- 
tex penetrates the cloud. Additional vortices then occur 
with increasing £1. We further analyzed how the vor- 
tices are affected by a tilt of the dipole axis. We find 
that although quantitatively the internal structure of the 
quantum state is in fact largely unchanged, the proba- 
blity density of the two-vortex state is sensitive to the 
tilt angle, very clearly mapping out the internal symme- 
try of the quantum state. From a theoretical perspective, 
this yields a fortunate situation as it allows the broken 
symmetry of the state to be analyzed directly via e.g. 
the density and current. For a rotationally symmetric 
Hamiltonian, this is typically not possible. In experi- 
ments, possibly any disorder or imperfection of the trap 
may already break the symmetry for aligned dipoles at 
(9 = 90.0°. 

In conclusion, our study suggests that there may be a 
possibility to experimentally observe vortices in a rotat- 
ing quantum system without superfluidity, composed of 
spin-polarized fermions with repulsive dipole-dipole in- 
teractions. 
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